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Abstract 

Using the high statistics r data sample recorded in the Belle experiment at KEKB, we have greatly 
improved the precision of the branching fraction B(t~ -> K *°(892)if^ T ) = (1.56 ± 0.02 ± 0.09) x 10" 3 , 
while the mass and width of the K*°(892) meson are measured to be (895.10 ± 0.27 ± 0.31) MeV/c 2 and 
(47.23 ± 0.49 ± 0.79) MeV/c 2 , respectively, with better accuracy than the PDG world average values. The 
first measurement of the decay r~ — > K* (892)K~tt°v t is also reported with a branching fraction of B(t~ — ► 
K*°(892)K-n iv T ) = (2.39 ± 0.46 ± 0.26) x 10~ 5 . 



INTRODUCTION 

In addition to high statistics samples of B 
mesons, B factories also provide large number 
of r-leptons, which can be used for high preci- 
sion measurements. We report here a study of 
the decay r~ -» K*°(892)K- v T , in which its 
branching fraction, as well as the mass, M K *o, 
and width, T^-.o, of the if* (892) resonance 1 
are determined. 

The measurement was performed with the 
Belle detector at the KEKB asymmetric- 
energy e + e~ (3.5 on 8 GeV) collider [1], using 
the reaction e + e~ — > t + t~. In this analysis 
we use a 544.9 fb _1 data sample recorded on 
and near the T(4S) resonance, corresponding 
to 5.0 x 10 8 produced r + r _ pairs. 

The Belle detector is a large-solid-angle 
magnetic spectrometer that consists of a sil- 
icon vertex detector (SVD), a 50-layer cen- 
tral drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintilla- 
tion counters (TOF), and an electromagnetic 
calorimeter (ECL) comprised of CsI(Tl) crys- 
tals located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An 



1 We hereafter simply denote the if* (892) meson as 
K*°. 



iron flux-return located outside the coil is in- 
strumented to detect K^ L mesons and to iden- 
tify muons (KLM). The detector is described 
in detail elsewhere [2]. 

Particle identification (PID) plays an im- 
portant role in this experiment and is based 
on a likelihood ratio V x = L x /^ y L y for a 
charged particle x = fx, e, K, or tt, where 
the sum runs over the relevant particles. L x is 
a likelihood based on the energy deposit and 
shower shape in the ECL, the momentum and 
dE/dX measured by the CDC, the particle 
range in the KLM, the light yield in the ACC, 
and particle's time-of-flight from the TOF. 

EVENT SELECTION 

We look for r~ -► K*°(892)K ~v T candi- 
dates with the following signature: 

Wal ~+ ^*°(892)+iv- + (missing) 

r tl g W e ) + + n(< 1)7 + (missing). 

An event should contain 4 charged tracks with 
zero net-charge. The following basic criteria 
are imposed; the total energy of the tracks 
and photons in the center-of-mass (CM) frame 
should be less than 11 GeV; the missing mo- 
mentum should be greater than 0.1 GeV/c and 
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lie within the detector acceptance —0.866 < 
cos 9 < 0.956. The event is subdivided into 
two hemispheres according to the thrust axis 
in the CM frame, the signal and tag sides. The 
signal side contains three prongs and no addi- 
tional photons while the one-prong side may 
include one extra photon. 

On the tag side, V^u > 0.1 is required 
for the track, while on the signal side the re- 
quirements V K > 0.8 (T> n > 0.8) with V e < 
0.9 are imposed for both kaons with opposite 
charges (for the K*° daughter pion) within 
cos 9 > —0.6. Here 9 is the polar angle with 
respect to the z axis, which is anti-parallel 
to the e + beam direction. The three sig- 
nal candidate tracks must have charge assign- 
ments that satisfy strangeness conservation i.e. 

In order to reduce qq background, the mo- 
menta of the K*° and K*°K~ system are 
required to satisfy the requirements p K ^> > 
1.5 GeV/c and p™ K > 3.5 GeV/c, respec- 
tively; the opening angles between the K*° 
and K~ and between the thrust axis and miss- 
ing momentum should satisfy the requirements 



cos ( 



K > °- 92 aIld COS Cust-miss < -0-6, 



respectively; the invariant masses of the parti- 
cles on both the tag and signal sides should be 
in the ranges M tag < 1.8 GeV/c 2 (~ m T ) and 
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Fig. 1: K +7r~ mass distribution after the selection. 
Data are indicated by points with error bars. The 
solid curve is the best fit result and the dashed 
curve is the incoherent BG N lnco (Mk^), see the 
text for details. 



M signal < 1.8 GeV/c 2 . 

After applying all the selections, 5.1 x 10 4 
events remain, and yield the K + ir~ mass 
(Mk-k) distribution, shown in Fig. 1. A clear 
K*° peak is seen. 

The detection efficiency (e) for 
t~ — ► K*°K~v T is evaluated using the 
KKMC/PYTHIA Monte Carlo (MC) pro- 
gram [3], where the V-A interaction is 
assumed for the weak vertices and all the 
final state particles are produced according 
to pure phase space. For K*° — > K + ir~, 
a spin-dependent Breit-Wigner (BW) func- 
tion is used. The resulting efficiency, e, is 
2.15% which includes the branching fraction 
B(K*° -► K + ir~) = 2/3 as well as the tagging 
efficiency. 



BACKGROUNDS 

The dominant background (BG) is from 
other generic decay modes of r lepton pairs. 
The largest component, composing ~80% of 
the BG, arises from t~ — > it + it~ K~nir°i> T 
and t~ — > 7r + 7r _ 7r _ n7r°i/ r with n =0 and 
1, through 7r + — > K + misidentification, and 
an undetected 7r°. We prepare such enriched 
BG samples from data by replacing the K + 
selection criteria Vk > 0.8 by Vk < 0.2 
(K + — ► 7r + ), and estimate the contamination 
in the K*°K~u T samples, taking into account 
the PID fake rate. This estimate has a 1.1% 
systematic uncertainty, mostly due to errors in 
the PID fake rate. 

Another ~20% of BG is attributed to r~ -> 
<l>K-v T , K + tt-K-tt°u t (excluding K*°), 
(/yir~v T , and K + ir~ K~v T (ex. K*°), and is es- 
timated by MC simulation. The first and sec- 
ond background components arise from mis- 
PID (K~ — > 7T _ ) and missing 7r° detection, re- 
spectively, at rates of (1.4 ± 0.2)% and (2.1 ± 
0.7)%. Their Mk-k spectra are evaluated by 
MC, based on their branching fractions, re- 
ported in Refs. [4] and [5], with uncertainties 
of 0.07% and 0.6% of the K*° yields. 

The third background component has 
the same final state, K + ir~ K~v T , as 
that of K*°K~v T , and its Mk-k spectrum 
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(N^ u (M K n)) is obtained in [6]. The fourth 
BG is the non-resonant (NR) component and 
can interfere with the signal. Details are 
discussed in the next section. 

On the other hand, i^*°-resonant BG comes 
from as yet unmeasured r~ — > K*°K~tt°u t 
decay and qq processes. In order to mea- 
sure B(t~ —> K* K~ir°v T ), we apply the 
same criteria as for the K*®K~v T selection, 
but with additional requirements of n 7 = 2 
and effective mass (M 77 ) in the range, 0.1178 
GeV/c 2 < M 77 < 0.1502 GeV/c 2 . Figure 2 
shows the resulting Mk-k distribution, where 
the K*° peak includes two contributions: one 
is the K*^K~y®v r signal, and the other comes 
from K*°K~u T events with a ir° formed by 
spurious photons, caused by a splitoff of a 
hadronic shower in the ECL calorimeter. The 
latter contribution is calculated by MC, us- 
ing B(t~ -> K*°K~u T ) from the PDG [5] 
and our result. The non-K*° BG in Fig. 2 
mostly arises from t + t~ pair processes, while 
the qq contribution is negligibly small. The 
M KlT distribution is fitted with a K*° BW 
function plus a non-K*° BG represented by 
a Landau function, with four free parameters: 
the total number of K*°K events, and 

three parameters of the Landau function. The 
best fit yields the number of signal events as 



N 



k*°k- 



129.2 ±25.1, with a detection ef- 

*0 j 



ficiency e = 0.54%, fixing the K K~v T peak- 
ing BG to N K *o K - u = 113. 7± 6.8 events. As a 
result, we obtain the first measurement of the 
branching fraction, 

K*°K-tt°ij t ) = 

(2.39 ± 0.46 ± 0.26) x 10~ 5 , (1) 



with a systematic uncertainty of 11.0%. Ta- 
ble II lists the sources of individual uncertain- 
ties. The largest contribution, 9.0%, arises 
from an uncertainty on the peaking BG esti- 
mate, especially, on r~ — > K*°K~v T . This 
K*°K~7T v T contamination in the K*°K~v T 
signal in the Mkk distribution is only 0.4%. 

The other possible i^*°-resonant BG, the 
qq contribution, mostly from uds, is examined 
in a comparison between data and MC, by ap- 
plying qq enriching selection criteria. The es- 
timated contribution in the K*° mass region is 



> 
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Fig. 2: K + ir mass distribution for K*°K 7r°t/ r 
candidates. Data are plotted by crosses, and the 
best-fit result is indicated by the solid curve, the 
K *°K~v T contamination and the non-K *° BG are 
shown by the dotted curve and dash-dotted one, 
respectively. See the text for details. 



46 ± 17 events: It gives an uncertainty of only 
0.14% to the estimation of the signal events for 
the K*°K~u T analysis. 



K*°K-v T ) 



The Mk-k distribution in Fig. 1 is fitted 
with the following formula: 

N(Mkk) = \oAbw(M k -k) + PAn^MkkY^ 

+N^~ u (M K k) + N incoh (M K K) (2) 

where Abw{Mk-k) is the Breit-Wigner (BW) 
function for the decay of a J p = 1~ — ► 0~0~ 
state, expressed as 



A bw (Mkk) 



M T 



r(M 



K-k = — 



(M 2 -M£j-iM r' 
\ 3 / Mp \ D L (q Q r) . 
qj \Mk-k) D L (qr)~ 



(3) 



Mo and To are the M K *o and r^*o, respec- 
tively, q is the momentum in the K + ir~ cen- 
ter of mass, and Di{qr) = 1/(1 + r 2 q 2 ) is 
the barrier factor with the so-called damp- 
ing factor r. Since the damping factors for 
if* (892) — > K + tt~ decay were obtained by 
LASS [7] in the K~p — ► Kirn reaction, r = 
3.4 ± 0.6 ± 0.3 (GeV/c) -1 , and by FOCUS [8] 
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Table I: Results of fits over the mass region of 0.645 — l.f47 (GeV/c 2 ) . 'standard' means the fit with 
five free parameters (a, f3, M^*o, rV»o and (ft) with r — 3.53, and with r = 0. The interference between 
Abw and Anr is excluded in the 'no interference' case, and then (ft is removed from the fit. In the 
rightmost column, Anr. is omitted in the fit. To obtain branching fractions, we use detection efficiencies, 
as ck*k — (2.15 ± 0.01)% and tKKix — (2.60 ± 0.01)%, determined from MC simulations based on pure 
hadronic phase space with a V — A weak interaction. In B{t~ — > K*°K~u T ), B(K*° — ► K + tt~) = 2/3 is 
included. 





standard 


no interference 


no A NR 




r = 3.53 r = 0.0 


r = 3.53 


r = 3.53 


M K *o (MeV/c 2 ) 
T K ,o (MeV/c 2 ) 
<M°) 

B(K*°K-is) (xlO- 3 ) 
B{K + t:-K-v) (xlO" 5 ) 


895.25 ± 0.27 896.58 ±0.28 
47.70 ±0.49 48.19 ±0.51 
63.46 ± 2.05 37.98 ±2.90 
1.56 ± 0.02 1.55 ±0.02 
5.76 ± 0.59 4.62 ±0.53 


896.72 ±0.19 
47.45 ± 0.48 

1.80 ±0.01 
4.82 ±0.54 


896.85 ±0.19 
49.21 ±0.46 

1.84 ±0.01 


X 2 /ndf 


83.29/80 84.64/80 


146.3/81 


226.2/82 



in D+ -» K-ir+fi+v decay, r = 3.96±0.54lJJJ 
(GeV/c)" 1 , we take the average of these two 
values, r = 3.53 ± 0.59 (GeV/c)- 1 in the 
fits. The term Anr(Mk-k) is the non-resonant 
— > K + ir~ K~v T scalar amplitude, and can 
be expressed as 

A NR = (^f 1 ) sm5 LASS e iSLASS (4) 

where the phase-shift of 5 lass is parametrized 
as 

r I bq 

cot d LA ss = f- — (5) 

aq 2 

following the LASS [7] and FOCUS [8] analy- 
ses with a = 4.03±1.72±0.06 (GeV)" 1 and b = 
1.29±0.63±0.67 (GeV)" 1 , obtained by LASS. 
The phase 5 lass has a weak g-dependence and 
is ~ 45° in our momentum range. Since the 
final hadronic system contains an additional 
K~ , final state interactions are possible be- 
tween K*° and non-resonant K + ir~ systems 
and K~, and a relative phase between the 
Abw and A^r amplitudes is introduced in the 
formula. Here, N^ w u (Mk-k) is the contribu- 
tion from <fiir~v T , and 7V mcoh (Mx7r) is the sum 
of the incoherent continuum and peaking BG's 
discussed in the previous section. These cor- 
respond to 795 ± 65 and 10, 073 ± 994 events, 
respectively. The free parameters in the x 2_ fits 
are a, (3, M K *o, T K *o, and 0. 



The mass dependence of the resolution is 
taken into account in the fit. The resolu- 
tion function can be approximated as a sum 
of three Gaussian functions with their relative 
fractions and standard deviations cr's: about 
52% with 1.7 MeV/c 2 , 38% with 3.4 MeV/c 2 
and 10% with 9.2 MeV/c 2 at M Kw ~ 890 
MeV/c 2 ; a's vary by 3.7 x 10" 3 , 5.6xl0~ 3 and 
3.3 x 10^ 3 per MeV/c 2 with M Kn , respectively; 
the third Gaussian is shifted to the higher mass 
side by ~ 0.13 MeV/c 2 relative to the others. 

The x 2 -fits are performed in the mass range 
Mktt = 0.645 - 1.147 (GeV/c 2 ). The re- 
sults of fits are listed in Table I, and the 
best fit is shown in Figs. 1 and 3. The 
branching fractions obtained are B{t~ — > 
K*°K~u T ) = (1.56±0.02) x 10~ 3 and B{t~ -» 
K+7r- J PC-^) non _ reS onant = (5.76±0.59) x 10- 5 . 
The non-resonant K + tt~ K~u t contribution is 
(5.58 ±0.57)%, which agrees well with the FO- 
CUS result of (5.30 ± 0.74±g;|§)%. 

Since <p + 5 L ass - 65° + 45° ~ 110°, the 
interference term 2oi(3AbwA nr has a Mk-k 
distribution similar to a BW shape, as seen 
in Fig. 3. The role of the damping factor in 
the Mk-k distribution in this reaction resem- 
bles that of <p in the interference term on the 
distribution. This can be seen in the = 0° 
case, where the parameters resulting from the 
fit, except for (ft, do not change much com- 
pared to the r = 3.53 case, while eft changes 
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M Kn (GeV/c 2 ) 



Fig. 3: Result of the best-fit, plotted with a loga- 
rithmic scale. The spectra of all terms, the inter- 
ference term and the non-resonant term are plot- 
ted by the solid, dashed and dotted curves, respec- 
tively. The negative interference term is reversed 
and plotted by the dash-dotted curve. See Fig. 1 
for the same plot, but with a linear scale. 

from (37.98 ± 2.90)° to (63.46 ± 2.05)° with a 
X 2 difference of A% 2 = 1.35. 

When the interference contribution is not 
taken into account, B(r~ — > K*°K~u T ) be- 
comes 15% larger (see the fourth column in 
Table I); when the NR contribution is totally 
ignored, £>(t~ — > K*°K~v T ) becomes 2.2% 
larger than in the case where the interference 
contribution is not taken into account (see the 
rightmost column in Table I). 

The systematic uncertainties are listed in 
Table II. The total uncertainty is 5.5%, 
where the largest errors are attributed to the 
track finding efficiency, and PID's on the lep- 
ton and kaon, as 3.3%, 2.9% and 2.6%, re- 
spectively. The r parameter is varied by 
±lcr to evaluate its effect on the branching 
fractions: the resulting change is ±0.3% for 
B(t~ -> K*°K-u T ) and ±1.6% for B{t~ -> 
K + tt~ K~v T ). These systematic uncertainties 
are much smaller than the statistical error of 
the fit. 

Adding all systematic errors in quadrature, 
the branching fraction obtained is 

B(t~ -> K*°K-u T ) = 

(1.56 ± 0.02 ± 0.09) x 10" 3 . (6) 

The branching fraction of non-resonant 



r — > K + tt K v T is also obtained, 

B(t > K TT K ZA-)non— resonant — 

(5.76 ± 0.59 ± 2.04) x 10" 5 . (7) 

This is the first measurement of non-resonant 
r~ — > K + it~ K~v r decay, where the large sys- 
tematic error mostly comes from the uncer- 
tainty in estimating the 7r + 7r~7r~z/ r BG con- 
tamination through mis-PID. 

M K ,o AND T K ,o 

The absolute mass scale for the Belle detec- 
tor is confirmed by reconstructing the K$ and 
<fi masses through their 7r~7r + and K~ K + de- 
cay modes, respectively, and comparing them 
to their world average values. Such studies 
provide mass differences from the PDG values 
[5] of -0.045 ±0.027 MeV/c 2 and -0.14±0.20 
MeV/c 2 for Mk s and respectively, or 

from the most precise measurements [9] and 
[10], -0.022 ±0.035 MeV/c 2 and -0.12 ±0.20 
MeV/c 2 . Therefore, no correction is applied 
to the absolute mass scale. 

On the other hand, MC simulation results 
in shifts of AM K * = 0.15 ± 0.04 MeV/c 2 and 
Ar x ,o = 0.47±0.22 MeV/c 2 , due to the effects 
of detection efficiency and event selections, and 
therefore we correct the raw M^*o and T^*o 
values by these amounts, respectively. 

Table II: Systematic errors for B(t~ -» K*°K~v T ) 
(left column) and B(t~ -> K*°K--k°v t ) (right col- 
umn) in %. 





K*°K- 


v K*°K-ir a v 


Luminosity 


1.4 


1.4 


cr(e + e _ — > tt) 


0.3 


0.3 


Tracking efficiency 


3.3 


3.3 


Trigger efficiency 


0.7 


0.1 


Lepton-ID 


2.9 


2.9 


Kaon-ID /fake 


2.6 


3.8 


MC statistics 


0.3 


0.5 


7T° efficiency 




1.7 


BG estimate 


1.3 


9.0 


Total 


5.5 


11.0 
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Table III: Comparison of our results with other experiments 





M K »o (892 ) (MeV/c 2 ) 


r*-.o (892) (MeV/c 2 ) 


Belle 
PDG [5] 
FOCUS [8] 
LASS [7] 


895.10 ±0.27 ±0.31 
896.00 ±0.25 
895.41 ± 0.32±g;i 
895.9 ±0.5 ±0.2 


47.23 ±0.49 ±0.79 
50.3 ±0.6 
47.79 ± 0.86±^§ 
50.8 ±0.8 ±0.9 




M K »± (S92) (MeV/c 2 ) 


IV± ( 892) (MeV/c 2 ) 


Belle [11] 
PDG [5] 


895.47 ±0.20 ±0.74 
891.66 ±0.26 


46.2 ± 0.6 ± 1.2 
50.8 ±0.9 



To evaluate uncertainties in the BG esti- 
mates on M K *o and T K *o, we vary the mis-PID 
probabilities by ±1<7 in the fit. The result- 
ing variation is 5M K *a = ±0.23MeV/c 2 and 
5Fk*o = ±0.75 MeV/c 2 . The uncertainty in 
r is included by varying it by ±lcr in the fit, 
which gives changes of M K *o = ±0.22 MeV/c 2 
and T K *o = ±0.09 MeV/c 2 . 

Including the above errors in quadrature, 
the mass and width are 

M K *o = 895.10 ±0.27 ±0.31 (MeV/c 2 ),(8) 
T K *o = 47.23 ± 0.49 ± 0.79 (MeV/c 2 ), (9) 

where the first and second errors are statistical 
and systematic ones, respectively. 

CONCLUSIONS 

With the use of high-statistics data sam- 
ples of 544.9 fb , collected in the Belle 
experiment, we have measured B(t~ — ► 
K*°K~v T ) with the highest precision so far 
attained, and obtained the first measure- 
ments of B(t~ -» K*°K-ir°v T ) and B(t~ 
K + Tr~K~i/ T ) uon - Te80na , at , respectively. 

The K*° mass and width are also mea- 
sured, respectively, with better or same accu- 
racy than those of the PDG world average val- 
ues, or the most precise previous experiments, 
as listed in Table III. 

Figure 4 compares various measurements: 
our results differ from the PDG world average 
values, but are quite consistent with the re- 
sults of the FOCUS experiment for M K *o and 
r^-*o. The difference between our measure- 
ment of B{j~ — > K*° K~v T ) and other results 



(a) Belle i-w-i 

I m i PDG 

ALEPH I • 1 

ARGUS ' ' * ' ' 

I 1 I I 

1.0 1.5 2.0 2.5 XlO 3 

Br(T->K*°Kv) 

(b) Be , le H • H 

I — • — I PDO 

FOCUS I I » — l-H 

LASS » • H 

892 894 896 898 

M K , (MeV/c 2 ) 

(c) i.i Belle 

I — • — I PDG 

FOCUS I 1 • 1 1 

LASS I — I • 1 — I 

I I I I I 

42 44 46 48 50 52 

r K . (Mev/c 2 ) 

Fig. 4: Comparison of various measurements: (a) 
B(t~ -> K*°K-v T ), (b) M K ,o and (c) r K .o. In 
(a), the PDG value is obtained from ARGUS [12] 
and ALEPH [13]. In (b) and (c), FOCUS and 
LASS results are included in the PDG value. 

may be due to the inclusion of the interference 
effect between the decays that have the same 
final state K + n~ K~v T . 
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